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ABSTRACT 

A method  to  simulate  microwave  heating  of  ceramics  which  has  a temperature 
dependent  dielectric  property  is  developed  here.  In  this  simulation  study,  the 
impedance  method  is  used  to  find  the  microwave  energy  absorbed  by  ceramics  and 
a non-linear  finite  element  method  is  used  to  determine  the  dynamic  temperature 
profile  in  the  ceramics  during  microwave  heating.  Using  the  developed  method,  the 
thermal  runaway  phenomenon  in  the  microwave  heating  of  ceramics  is  successfully 
simulated.With  detailed  analysis  of  the  microwave  energy  absorption  pattern  in  the 
ceramics,  the  effects  of  dielectric  properties  on  microwave  energy  absorption  by 
ceramics  are  discussed.  The  causes  of  non-uniform  heating  with  microwave 
energy  that  has  been  observed  in  our  laboratory  are  also  investigated. 

INTRODUCTION 

The  use  of  microwave  energy  is  a new  and  exciting  approach  in  ceramic 
processing.  It  has  already  been  used  in  sintering,  joining  and  melting  of 
ceramics[l].  Since  microwave  heating  is  a volumetric  process,  it  could  provide 
uniform  heating  so  that  the  temperature  gradient  vhich  is  observed  in  conventional 
rapid  heating  methods  can  be  avoided.  Rapid  and  uniform  heating  are  important  in 
the  joining  and  sintering  of  ceramics.  On  the  contrary,  non-uniform  heating  is  often 
observed  in  our  laboratory  with  microwave  sintering  or  joining.  Therefore,  it  is  of 
practical  interest  to  simulate  the  phenomenon  of  microwave  heating  for  better 
control  and  more  efficient  use.  In  spite  of  the  significance  of  the  problem,  there  is 
no  comprehensive  analysis  available  which  would  describe  the  behavior  of  ceramic 
materials  exposed  to  electromagnetic  radiation.  Research  by  Iskander  [2]  and 
Watters  et  al.  [3]  has  revealed  some  of  the  mechanisms  of  microwave  heating  of 
ceramics.  However,  the  simulation  of  microwave  heating  of  ceramics  with  a 
temperature  dependent  dielectric  property  is  still  lacking.  In  this  paper,  a method  of 
simulating  microwave  heating  of  ceramics  with  temperature  dependent  dielectric 
properties  is  developed  here.  The  impedance  method  is  used  to  find  the  microwave 


energy  absorbed  by  ceramics.  A non-linear  finite  element  method  is  developed  to 
determine  the  dynamic  temperature  profile  in  the  ceramics  during  microwave 
heating.  Using  this  'Jcveloped  method,  the  thermal  runaway  phenomenon  in  the 
microwave  heating  o.  ceramics  is  successfully  simulated. With  detailed  analysis  of 
the  microwave  energy  absorption  pattern  in  the  ceramics,  the  effects  of  dielectric 
properties  on  microwave  energy  absorption  by  ceramics  are  discussed.  The 
causes  of  non-uniform  heating  using  microwave  energy  that  has  been  observed  in 
our  laboratory  are  flso  investigated.  In  doing  so,  a better  understanding  of 
microwave  heating  of  ceramics  is  realized. 

THEORY 

In  order  to  simulate  microwave  heating  of  ceramics,  it  is  necessary  to  find  the 
electric  and  magnetic  field  strength  inside  ceramics  and  the  absorbed  microwave 
energy.  Electric  and  magnetic  fields  are  linked  by  Maxwell’s  equations,  a group  of 
linear  differential  equations.  Assuming  an  e_ia)t  harmonic  time  dependence, 
Maxwell’s  equations  can  be  expressed  as  follows. 


V • ) = pe 

(1) 

V.  (HoUrH^O 

(2) 

'*  x E = j coPoPfH 

(3) 

Vxtl  = aE-j 

(4) 

where  e0  and  |i0  are  the  permitivity  and  permeability  in  the  vacuum,  e,.  and  (ij  are 

the  relative  permitivity  and  permeability  of  the  material , a is  the  conductivity  of  the 
material.  E and  H are  the  electric  and  magnetic  field  strength,  respectively.  The 
propagation  of  energy  in  the  electromagnetic  field  can  be  deduced  from  this 
equation  system  and  leads  to  Poynting’s  theorem 


P = - J[ExH*  ]dS  (5) 

s 

which  states  that  the  mean  energy,  P,  flowing  into  a surface,  S,  depends  on  the 
amplitude,  distribution  and  prevailing  phase  of  the  electric  and  magnetic  field. 

By  using  Gauss’  law,  equation  5 can  be  converted  into  the  volume  integral  which 
can  then  be  resolved  into  three  single  integrals 
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(6) 


(e  • E * )dv 


The  first  two  integrals  take  account  of  the  magnetic  and  electric  fields  respectively 
while  the  third  represents  the  necessary  dissipation  in  the  dielectric  in  a general 
form.  Therefore,  the  energy  converted  into  heat  by  the  alternating  field  is 


P=cae0£j|(E*E*)dv  (W)  (7) 

v 

It  increases  according  to  the  frequency , the  square  of  the  electric  field  strength  and 
the  imaginary  part  of  the  dielectric  constant  Once  the  profile  of  e as  function  of 
temperature  and  the  electric  field  strength  in  the  homogeneous  body  are  known,  it 
may  be  possible  to  study  the  thermal  runaway  conditions  through  the 
source-incorporated  heat-diffusion  equation.  The  diffusion  of  thermal  energy  in  a 
homogeneous  bound  volume  V is  determined  by  the  partial  differential  equation 

PCpff-V'l'P  (8) 

where  p,  Cp  and  K|,  are  the  mass  density , specific  heat  and  thermal  conductivity  of 
the  material,  respectively,  p is  the  microwave  energy  density  absorbed  by  the 
material.  At  the  boundary  of  the  volume  V,  the  boundary  condition 


Khn.VT  = h(T-T0)  + ear(T4-T04]  (9) 


must  be  satisfied , in  which  h is  the  he*:  convection  coefficient,  e and  ar  are  the 
emissivity  of  the  material  ana  Stefui-Boitzmann  constant.  T0  is  the  ambient 
temperature.  The  initial  condition  is 


T(r , 0)  = Ti  (10) 

This  heat  diffusion  equation  is  analogous  to  the  forced  Fisher  equation 

Tt  - T + G ( T ) (11) 
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which  is  known  to  have  chaotic  behavior  for  specific  initial  and  boundary 
conditions  as  investigated  by  Fisher  [ 4 ] in  1937  and  Rothe  [5]  in  1981. 


DESCRIPTION  OF  THE  MODEL 

For  simplicity,  a ceramic  slab  with  finite  thickness  under  plane  wave  radiation 
is  considered  here  as  depicted  in  figure  1.  The  incident  electric  field  is  a 
monochromatic  plane  wave  propagatir  g in  tee  z-direction  and  is  polarized  along  the 
x-axis.  To  account  for  material  non-lineanty  during  microwave  heating,  the  slab  is 
further  divided  into  layers  so  mat  the  finite  element  method  can  be  used  accurately. 
It  is  assumed  that  each  element  has  the  same  material  properties  during  the 
microwave  heating  process  at  each  temperature  step.  Since  the  ceramic  slab  is 
assumed  to  be  very  large,  the  problem  becomes  one-dimensional.  Li  the  following 
discussion,  layers  with  smaller  thicknesses  will  be  cons’dered  as  ca'ferenc  meoia 
since  they  may  have  different  material  properties  suDn  as  dielectric  constant  and  loss 
factor  which  are  functions  of  temperature  during  the  inicrowavt  heating  process. 


Figure  1.  Cerarmc  Slab. 

To  simulate  microwave  heating  of  ceramics,  the  microwave  energy  absorbed  by 
the  ceramic  slab  must  first  be  calculated.  Therefore,  the  electric  and  magnetic  fields 
in  the  ceramic  slab  have  to  be  determined  first.  The  electromagnetic  field  inside  a 
dielectric  body  of  arbitrary  shape  is  difficult  to  determine.  For  the  model  considered 
here,  ar.  impedance  method  [6]  is  applied  and  proven  to  be  effective  to  account  for 
material  non-lineanty. 

The  application  of  the  impedance  method  can  be  described  a;  folio  vs.  When  an 
electric  field  is  incident  on  the  ceramic  slab,  which  is  now  assumed  to  be  composed 
of  several  different  layer,  multiple  leflection  will  occur  in  the  different  !ayers 
leading  to  positive  and  negative  traveling  waves.  A generalized  reflection  coefficient 
is  defined  for  any  layer  as  the  ratio  of  the  n.cident  and  reflected  field.  Hence,  the 
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totai  field  can  be  renresented  by  tne  positive  traveling  wave  and  tne  generalized 
'eriecuon  coerficiem . in  the  first  layer,  the  incident  field  is  known  and  the  reflected 
field  is  unknown.  For  the  Nth  layer,  tne  generalized  reflection  coefficient  is  zero 
because  there  is  no  negative  traveling  wave  The  total  field  impedance,  which  is 


mpiex  and  position  dependent,  is  defined  a • the  ratio  of  electric  field  over 
ignetic  field.  For  tne  Nth  iaver.  tne  total  fief  ' J T1  ~ 


impedance  will  oe  equal  to  the  field 
impedance  tor  met  layer.  At  an  m terrace,  the  totai  field  impedances  for  adjacent 
'avers  are  eaual  due  to  the  field  boundary  condition  at  that  interface.  With  this 
-ouncary  conditior.  the  totai  impedance  at  different  interfaces  can  be  found.  Hence, 
re  re:iected  field  m the  first  layer  is  found.  By  using  the  boundary  condition  that 
re  electric  field  is  continuous  at  each  interface,  the  electromagnetic  field  strei.stn  in 
_ach  layer  :s  readily  found.  The  resulting  field  can  men  be  used  to  calculate  the 
microwave  energy  aosorbed  by  that  layer. 

When  'he  microwave  energy  absorbed  bv  the  slab  is  known,  tne  heat  diffusion 


■non  can  oe  used  to  caicuia 


r ai  ma  JAm 


nperature  variation  with  time  and  position. 


Since  me  dielectric  constant  and  loss  factor  are  functions  of  temperature,  the 
microwave  energy  absorbed  by  the  ceramic  siab  is  aiso  a function  of  temperature. 
Hence,  the  heat  diffusion  equation  becomes  non-iinear.  A non-iinear  finite  element 
method  is  therefore  reeded  to  find  the  dynamic  temperaiuie  distribution  profile.  At 
the  surface  of  the  ceramic  siab.  radiation  link  elements  are  used  to  account  for 
radiation  heat  loss.  Conduction  loss  is  neglected  since  the  radiation  loss  is  me  pnme 
heat  loss  at  high  temperature.  The  detailed  imt  lamentation  of  the  non-iinear  analysis 
is  is  follows.  The  time  step  to  do  the  non-linear  analysis  is  designated  first,  the 
ter  merature  at  the  end  of  time  step  is  then  estimated,  the  material  properties  at  the 
middie  of  the  temperature  increment  are  used  for  each  media.  The  microwave 
energy  absorbed  by  each  media  with  different  dielectric  properties  is  calculated 
according  o the  technique  described  above.  The  temperature  distnbuticn  is  then 
obtained  by  using  power  absorption  data.  The  computed  resuits  will  be  compared 
with  the  p-or  estimated  temperature.  Such  an  iteration  procedure  will  continue  until 
■he  difference  between  the  estimated  and  calculated  temperature  reaches  a prescribed 
value. 

RESULTS  AND  DISCUSSION 
Microwave  Energy  Absorp..on  by  Ceramics 

By  using  the  echnique  stated  above,  the  effects  of  the  dielectric  constant  and  the 
loss  factoi  on  the  power  absorption  by  ceramics  ure  considered.  Figure  2 gives  the 
comparison  of  rower  absorption  by  ceramic  slabs  with  the  s^me  dielectric  constant 
and  different  loss  factors.  Hie  increase  in  ioss  factor  will  dramatically  increase  the 
ability  of  the  ceramic  siab  to  absorb  microwave  energy,  figure  5 shows  that  with 
the  increase  of  bom  dielectric  constant  and  loss  factor,  the  uniformity  and  ability  of 
power  absorption  by  ceramics  are  aiso  increased. 

Simulating  Microwave  Heating  of  Ceramics 


: r.e  developed  r.on-imear  .unite  element  metnod  is  used 
tvnamic  temoeratuie  profile  01  a ceramic  slab  under  t 


h;s  case  tc  find  the 
■e  wave  radiation 


Micro1 


considering  changes  in  both  the  dielectric  constant  and  the  loss  factor  with 
temperature.  The  analysis  procedure  is  displayed  h figure  4.  The  data  of  dielectric 
constant  and  loss  factor  change  with  temperature  are  taken  from  Fukushima  et 
al.(7].  The  incident  microwa  ve  power  flux  is  30kw/m2.  The  microwave  frequency 
is  taken  to  be  6 GHz  to  be  consistent  with  the  dielectric  data. 
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Figure  2.  Microwave  Power  Absorption  for 
Slabs  with  Different  Loss  Factors. 


( Incident  Wavelength  ) 

Figure  3. Microwave  Power  Absorption  for  Slabs 
with  Different  Dielectric  Properties. 
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The  mass  density,  specific  heat  and  thennai  conductivity  are  taken  to  be  4g/cmA 
1.125  j/g  °C  and  iOw/cm2,  respectively.  Tie  thickness  of  the  plate  is  taken  to  be 
5.08  cm.  The  finite  element  analysis  routine  on  AN, SYS  is  implemented  in  the 
calculation.  The  calculation  is  done  on  a VAX- 1 1/780  computer.  The  CPU  time  is 
79  seconds. 


Figure  4.  Procedures  for  Simulating  Microwave  Heating  of  Ceramics. 

Figure  5 gives  the  temperature  variation  with  time  at  0.01m  inside  the  slab  from 
the  microwave  incident  plane.  7r ie  temperature  me ’esses  slowly  at  the  beginning 
and  rapidly  after  600°C. 

Figure  6 displays  the  temperature  profile  over  the  thicknest  of  the  slab.  An 
appreciable  temperature  gradient  is  observed.  Thi.  temperature  gradient  is  caused 
by  non-symmetric  microwave  radiation  of  the  ceramic  slab  and  rad.ation  heat  loss  at 
the  boundary,  which  subsequently  results  in  a non-uniform  power  absorption  by 
the  ceramic  slab.  If  symmetric  radiation  ,s  realized,  i.e.,  microwave  radiation  is 
from  both  sides  of  the  slab,  the  uneven  heating  will  result  from  boundary  radiation 
heat  loss  only  and  the  center  of  slab  will  have  the  highest  temperature.  Hence,  the 
radiation  heat  loss  at  the  bound-cv  is  the  main  contribution  to  the  non-uniform 
Meeting  wrih  microwaves  observec  in  our  laboratory  where  the  ceramic  sample  is 
melted  af  the  center  while  the  boundary  is  still  intact.  In  ordei  to  prevent  this  etfect. 
gooa  insulation  must  be  used  at  the  boundary.  Also,  in  practice,  the  ceramic  sample 
needs  to  b:  rotated  continuously  to  pre  rent  any  uneven  radiation.  Figure  7 shows 
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the  variation  of  total  microwave  energy  absorbed  by  the  ceramic  slab  verse  the 
temperature  variation  at  0.01m  inside  the  slab  from  microwave  incident  plane.  It 
indicates  that  as  the  ceramic  becomes  hot,  its  energy  absorption  ability  is  increased. 
Therefore,  thermal  runaway  is  realized  in  microwave  heating. 


Figure  5.  Dynamic  Temperature  Profile  O.Olm 
inside  the  Slab  from  the  Incident  Plane. 


Figure  6.  Temperature  Distribution  over  the  Thickness  of  the  Slab. 
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Figure  7.  Total  Microwave  Energy  Absorbed  by  slab  vs.  Temperature. 
CONCLUSION 

A method  of  modeling  microwave  heating  of  ceramics  is  developed  here.  The 
results  show  that  increasing  the  dielectric  constant  could  increase  microwave  power 
absorption  uniformity  while  increasing  the  loss  factor  could  increase  the  material's 
ability  to  absorb  microwave  energy.  It  is  found  that  non-uniform  heating  observed 
in  the  laboratory  can  be  caused  by  boundary  radiation  loss  and  non-uniform 
radiation  by  the  microwave  source.  Through  this  research,  it  is  observed  that  the 
dielectric  property  of  a material  at  elevated  temperature  has  a very  important  role  in 
designing  microwave  processing  technique.  In  microwave  sintering  of  ceramics, 
the  green  sample  changes  its  microstructure  during  microwave  heating.  Its  dielectric 
properties  change  with  not  only  temperature  but  also  microstructure.  Hence,  the 
characterization  of  the  dielectric  property  of  ceramics  during  microwave  sintering  is 
very  important . In  doing  so,  we  are  not  only  able  to  control  the  sintering  process 
but  also  able  to  understand  thoroughly  the  mechanism  of  microwave  sintering. 
Therefore,  we  need  to  develop  a method  to  dynamically  characterize  the  dielectric 
property  of  materials  with  temperatures  as  well  as  model  microwave  heating  for  the 
ceramic  samples  of  complicated  shapes. 
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